Lensed CMB simulation and parameter estimation 
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Modelling of the weak lensing of the CMB will be crucial to obtain correct cosmological parameter 
constraints from forthcoming precision CMB anisotropy observations. The lensing affects the power 
spectrum as well as inducing non-Gaussianities. We discuss the simulation of full sky CMB maps 
in the weak lensing approximation and describe a fast numerical code. The series expansion in the 
deflection angle cannot be used to simulate accurate CMB maps, so a pixel remapping must be 
used. For parameter estimation accounting for the change in the power spectrum but assuming 
Gaussianity is sufficient to obtain accurate results up to Planck sensitivity using current tools. A 
fuller analysis may be required to obtain accurate error estimates and for more sensitive observations. 
We demonstrate a simple full sky simulation and subsequent parameter estimation at Planck-like 
sensitivity. The lensed CMB simulation and parameter estimation codes are publicly available. 



I. INTRODUCTION 



The CMB temperature and polarization anisotropies 
are being measured with ever more precision. The statis- 
tics of the anisotropies can provide valuable limits on cos- 
mological parameters as well as constrain early universe 
physics. As we enter the era of precision measurement 
with signal-dominated observations out to small scales, 
the non-linear effects will become important. One of the 
most important of these on scales of most interest for pa- 
rameter estimation is that of weak lensing. Fortunately 
it can be modelled accurately using linear physics: the 
linear potentials along the line sight lensing the linear 
perturbations at the last scattering surface 0,H,|i|. Mod- 
elling of fully non-linear evolution is not required for the 
near future on scales with I < 2000, and non-linear cor- 
rections can be applied to the lensing potential if and 
when required. On smaller scales the situation becomes 
much more complicated anyway due to point sources, 
beam size and other non-linear effects. 

Lensing induces non-Gaussianities in the lensed CMB 
sky, and also changes the power spectra of the pertur- 
bations. Lensing will start to have an observable effect 
on the power spectrum very shortly, and hence needs to 
be taken into account to obtain correct parameter con- 
straints and error bars. For future observations, includ- 
ing the Planck 1 satellite and forthcoming ground based 
telescopes the effect is very important. In this paper 
we describe the simulation of lensed CMB maps (includ- 
ing the full non-Gaussian structure), but show that us- 
ing an accurate calculation of the lensed CMB power 
spectra |^| a naive parameter estimation (neglecting non- 
Gaussianities) works rather well up to Planck sensitivi- 
ties. Observations at higher sensitivities and resolutions 
will probably require a fuller analysis accounting for the 
full non-Gaussian distribution of the lensed sky, an im- 
portant problem that we do not tackle here. Our simu- 



lation code can be used for testing future methods, and 
the simple power spectrum parameter estimation method 
can act as a useful baseline for future improvements. 



II. WEAK LENSING OF THE CMB 

The small scale CMB anisotropy is dominated by the 
emission from the last scattering surface at redshift z ~ 
1000. Weak lensing of the CMB deflects photons coming 
from an original direction n' on the last scattering surface 
to an observed direction n on the sky today, so a lensed 
CMB field is given by 



X(n)=X(n') 



(1) 



where X is the unlensed field. The arcminute-scale dis- 
placement of the points is determined by the potential 
along the line of sight to the last scattering surface, con- 
veniently encapsulated into an integrated lensing poten- 
tial ip Q. The deflection vector is given by the gradi- 
ent of the lensing potential V?/'(n), where V is the co- 
variant derivative on the sphere. The vector n' is ob- 
tained from n by moving its end on the surface of a unit 
sphere a distance |V^(n)| along a geodesic in the di- 
rection of V?/>(n) [3, |5j. This is sometimes written as 
n' = n + V?/>(n). To our level of approximation |V^| is 
assumed to be constant between n and n', consistent with 
working out the lensing potential in the Born approxi- 
mation (evaluating the potential along the undeflected 
path). Lensing deflections are a few arcminutes, but are 
coherent over degree scales, so this is a good approxima- 
tion. 

The power spectrum for the lensing potential (jf (and 
the correlation to the temperature Cf T ) can be com- 
puted numerically in linear theory for a particular cos- 
mological model using CAMB 2 0, [a, • From the power 
spectrum simulated realizations can be made assuming 



*URL: http://cosmologist.info 
1 http : //sci . esa. int/pla nck| 
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FIG. 1: The effect of lensing on the CMB power spectra. The top plots show the fractional change in the temperature spectrum 
Cf T and the lensing-induced B-polarization spectrum Cp B . The bottom plots show the lensed (grey/red, less peaked) and 
unlensed (black) T-E cross-correlation Cj B and E-polarization Cf E power spectra. All results are for the fiducial model given 
in the text, and the lensed B-mode power spectrum shown is not very accurate due to the neglect of non-linear evolution in 
the lensing potential. 



Gaussianity of the primordial fields. Non-linear evolu- 
tion of the potential changed the power spectrum Cff 
on small scales and also makes the ip distribution some- 
what non-Gaussian. The power spectrum l 4 Cf^ peaks 
at Z ~ 60, however small scales where non- linear evolu- 
tion is important at late times only contribute to ip at 
I 3> 60. On these scales the contributions to the lens- 
ing potential come from a rather broad range of redshifts 
from 1 < z < 10, so even at I ~ 1000 the late time 
non-linear evolution does not radically change Cj . We 
therefore neglect the small effect of non-linear evolution 
here, though it can become important for high resolu- 
tion polarization B-mode experiments. Using CAMB with 
halofit 0, Q the non- linear evolution of if) can be esti- 
mated to change the lensed temperature power spectrum 
Cj T by about ~ 0.2% at I ~ 2000, though growing to 
1% or more on smaller scales. 

The significant effect of lensing on the CMB power 
spectra is shown in Fig. ^ where the lensed power spec- 



tra are computed accurately numerically using the corre- 
lation function method of Ref. . The lensing smoothes 
out features in the temperature and polarization power 
spectra, changing the C; peaks by up to 20% for the E- 
polarization on the scales of interest. Weak lensing does 
not change the total variance of the CMB anisotropics, 
with 

j2(2i+i)Q = j2( 2l + 1 ) c i' ( 2 ) 
i i 

and similarly for the polarization |^|. This encapsulates 
that fact that the observation in any fixed direction is 
just a displaced view of the last scattering surface, and 
hence is Gaussian and with the same variance as if there 
were no lensing. Lensing induces a non-Gaussian spatial 
correlation structure to the lensed CMB fields but does 
not alter the variance at a point. 

We now move on to discuss how to simulate full-sky 
lensed CMB maps. 
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FIG. 2: Power spectra from realizations of the 1st, 2nd, 3rd 
and 4th order terms in the lensing potential series expansion 
of the lensed temperature T; m compared to the full lensed Ci 
(top). The spectra may contain some pixelization error. 



to proceed is to re-map the anisotropy field by the deflec- 
tion vector V-0 as a function of position. The deflection 
field is a vector field on the sphere, and can be easily sim- 
ulated using spin-1 spherical harmonics (or, equivalently, 
vector spherical harmonics). Full details of the simula- 
tion process for the temperature and polarization on the 
sphere are given in the Appendix, along with details of 
a fast multi-processor implementation using a modified 
version of HEALPIX 3 1.2 |lCj . 

It is important to remember that the lensing deflects 
the physical anisotropy field, not the field after beam 
and pixel convolution. For this reason sky simulation 
and re-mapping have to be performed at high resolution 
regardless of how broad the observational beam is. The 
lensed field can then be convolved as required. 



III. PARAMETER ESTIMATION 



A. Series expansion 

A common procedure for working with CMB weak 
lensing it to perform a scries expansion in the deflection 
angle, so for the temperature 

T(n) = T(n + VV>) 

= t + v a rv a ^ + |v a vvVv a v fc r + ...(3) 

and the last line is evaluated at n. This expansion will 
be valid for scales much smaller than the deflection an- 
gle. On smaller scales the unlensed fields are deflected 
by a distance comparable to their wavelength, and the 
change in phase is not a small perturbation. Even on 
scales where the expansion is valid, the expansion only 
converges relatively slowly. The power spectra of the in- 
dividual terms in the expansion are shown in Fig.[5]up to 
fourth order. This clearly shows that a low order series 
expansion cannot be used for accurate map simulation 
due to the large variance of higher order terms on small 
scales. The lowest order series expansion may be a useful 
approximation for some applications (e.g. Ref. 9]), but 
can easily give results which are sufficiently inaccurate 
to be problematic. For example it was shown in Ref. 4] 
that the lensed CMB power spectrum computed from the 
lowest order series expansion of Ref. \J\ gives lensing cor- 
rections which are incorrect by an order unity factor on 
small scales. 



B. Re- mapping points 

Since the series expansion is not accurate enough for 
making simulated maps at high resolution, the best way 



With future precision observations we would ideally 
like to extract information about both the unlensed CMB 
(which has simple Gaussian statistical properties), and 
about the lensing potential (which contains additional 
information). Gravitational lensing of scalar E-mode 
polarization may ultimately dominate the B-mode po- 
larization signal from primordial gravitational waves, so 
an accurate treatment of the lensing will become essen- 
tial [Til IT2I IT3I ]. In the more immediate future grav- 
itational lensing will have a significant effect on the 
statistics of the observed CMB, and must be accounted 
for to obtain reliable parameter estimates. The non- 
Gaussianities induced by lensing must also be accounted 
for when attempting to assess the degree of primordial 
non-Gaussianity (an important probe of early universe 
physics) . 

The amount of information that can be learned about 
the lensing potential from the non-Gaussianities depends 
on the noise level [SJ llJ, [Tjj- For observations up to 
Planck sensitivity any reconstructed map of the lensing 
potential would be completely noise dominated, and the 
extra information this contains is rather limited. How- 
ever the effect of lensing on the power spectrum is sig- 
nificant and certainly cannot be ignored. In this paper 
we do not address the problem of handling the full like- 
lihood function, but merely show that using the lensed 
power spectrum and approximating the lensed field as 
Gaussian is sufficient to obtain good parameter estimates 
at Planck sensitivity. More general methods for handling 
the lensing likelihood function will be discussed in a fu- 
ture paper if they can be made to work. 



http : //www . eso . org/science/healpix/ 
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A. Gaussian Ci Likelihood function 

Spherical harmonic coefficients of a Gaussian temper- 
ature and polarization field on the sky can be used to 
define estimators of the covariance 



(4) 



where Wi m and Xi m are spherical harmonic coefficients 
of the temperature T, E-polarization E or B-polarization 
B. These are unbiased in that averaged over realizations 
{C™ X ) = C™ X . 

The assumed Gaussianity of Tz m , Ei m and Bi m gives 
the following full-sky likelihood function: 



-21ogP(Q|Q) = (21+1) {n 



log|C ; |} (5) 



(to within an irrelevant constant) where 







C BB 



(0) 



and C is the corresponding matrix of estimators. In the 
presence of instrumental noise the Ci and Ci should in- 
clude the noise variance. We have assumed a statistically 
parity invariant ensemble so that C BT — C BE = 0. 

For the lensed sky this is not the correct relation be- 
cause the lensed sky is not Gaussian if the lensing poten- 
tial is not fixed. However as discussed below, replacing 
Ci by Ci does not give significantly biased results for 
Planck, and is by far the simplest method of account- 
ing for CMB lensing in a parameter analysis: i.e. just 
pretend the lensed field is Gaussian, and use the theo- 
retical lensed CMB power spectrum. However at high 
resolutions and sensitivity this will not be correct, as 
the statistics of the Ci on small scales are governed by 
the same small number of large scale lensing modes [l6| . 
The naive approach could be improved by (for exam- 
ple) calculating an empirical Ci covariance from simula- 
tions and using this to make some partial correction for 
non-Gaussianity induced variations to the posterior dis- 
tribution. Non-Gaussian corrections should be small at 
Planck sensitivity, however the non-Gaussian corrections 
to futuristic signal dominated lensing B-mode observa- 
tions can be important ^tJ- We have checked that the 

lensed estimator Ci TT and Ci EE variances from simu- 
lations agree with the Gaussian results to within a few 
percent at I < 2000. 



B. Sampling from the posterior 



power spectrum with spectral index n s and amplitude 
A s (at wavenumber 0.05Mpc _1 ), baryon density Fl^h 2 , 
cold dark matter density £l c h 2 , optical depth r (reion- 
ization assumed sharp), and Hubble parameter today 
Ha = lOO/i km s -1 Mpc . We approximate the neutri- 
nos as massless and assume standard general relativity. 

We use the CosmoMC 4 jj] Markov Chain Monte 
Carlo (MCMC) package to sample from the posterior 
distribution of the parameters given the observed (simu- 
lated) data. To make the posterior parameter distribu- 
tions more Gaussian we use 9 r as a base parameter (with 
flat prior) instead of the Hubble parameter Ho. The de- 
rived parameter 9 r is defined as the (approximate) ratio 
of the sound horizon at last scattering to the angular di- 
ameter distance |l9j. a non-linear function of the other 
parameters that is very well constrained by the position 
of the acoustic peaks. We also transform to the am- 
plitude parameter ]xiA s (with a flat prior) that is con- 
strained very well in a linear combination with r because 
A s e~ 2r determines the small scale amplitude of the Cj. 
We then use the covariance matrix to transform to an 
uncorrelated set of parameters that the MCMC proposal 
density can use to explore the posterior distribution ef- 
ficiently. Given an approximate covariance matrix from 
previous runs chains converge in a few hours. For the 
first run the covariance can be learned dynamically as 
the chain evolves (discarding samples from the evolving 
part of the chain), or one can use the Hessian at the best 
fit point as a useful starting approximation. 

Since the computation time for generating the lensed 
Ci is dominated by the time to compute the transfer func- 
tions for the unlensed Ci and the lensing potential, pa- 
rameters like lnA s and n s remain 'fast' parameters [18j, 
in that changing them is quick as long as the other param- 
eters remain fixed. Thus methods to efficiently exploit 
the difference between 'fast' and 'slow' parameters [2£j 
can still be used to speed up MCMC runs even when 
CMB lensing is included. 

We neglect C BB since it is noise dominated even at 
Planck sensitivity and has almost no effect on parameter 
constraints, though including it in the Gaussian approx- 
imation is trivial if desired. 



C. WMAP data 

The importance of the lensing effect depends on the 
amplitude of the potential fluctuations along the line of 
sight. Larger amplitudes cause more lensing. By itself 
the first year WMAP data constrains the amplitude 
rather poorly due to a degeneracy with the optical depth 
and the absence of Cf E data to give a good upper limit. 
Large values of the optical depth r > 0.5 allowed by 
the data correspond to models with large amplitudes in 



We assume a simple flat adiabatic KCDM cosmolog- 
ical model with the following parameters to be deter- 
mined from the data: primordial curvature perturbation 
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which the lensing effect is quite significant, so strictly 
speaking the lensing should be accounted for in analyses 
of the first year WMAP data alone [HEl- In a totally 
free WMAP first year parameter analysis the lensing al- 
ready has a noticeable effect. However the region of pa- 
rameter space with r > 0.5 is almost certainly disallowed 
by numerous other sources of data and so should not be 
taken too seriously. 

The second year data WMAP should include C EE 
and a better measurement of the third peak, restricting 
to a much smaller range of r and flbh 2 , and so mak- 
ing the effect of lensing on the tails of the distribution 
less significant. As a simple toy model consider a full 
sky observation with isotropic Gaussian noise with vari- 
ance Nf T = N EE /A = N t BB /4: = 0.03/iA 2 and a sym- 
metric Gaussian beam of 13 arcminutes full- width half- 
maximum and neglect foregrounds. We find that neglect- 
ing the lensing effect on the power spectrum leads to a 
posterior mean of f^/i 2 about half a standard deviation 
lower than it should be for a typical realization. This is 
easily understood: lensing smoothes out the third peak 
in the temperature C/, making it appear lower relative 
to the first and second peaks by a couple of percent. A 
similar effect can be produced without lensing by low- 
ering the baryon density, so unless lensing is modelled 
consistently there is a danger of confusion. 

Future ground or balloon results that can resolve the 
higher temperature C; peaks will also be sensitive to lens- 
ing as the corrections become larger than 5% at I > 1500. 
We conclude that very soon it will be important to in- 
clude lensing in parameter analyses to obtain accurate 
results from CMB observations. 



D. Planck-like simulation 

As an example of a future observation where mod- 
elling the lensing will be crucial we now consider a sim- 
ple full-sky CMB observation simulation at Planck-like 
sensitivity. We compute theoretical CMB power spectra 
for the lensing potential and unlensed fields for a sim- 
ple fiducial model with n s = 0.99, A s — 2.5 x 10~ 9 , 
Q b h 2 = 0.22, Q c h 2 = 0.12, t = 0.15, h = 0.72, and sim- 
ulate full sky lensed maps with 12 x 2048 2 ~ 5 x 10 7 
pixels (generated by remapping a 12 x 8192 2 ~8x 10 8 
pixel unlensed sky as described in the appendix). As a 
simple toy model of an observation at optimistic Planck- 
like sensitivity we assume isotropic Gaussian noise with 
variance Nf T = N EE /4 = N BB /A = 2 x 10~ 4 /iA 2 and 
a symmetric Gaussian beam of 7 arcminutes full-width 
half-maximum and neglect foregrounds. 

We use simulated lensed Ci on scales with I < 2000. 
As discussed above we neglect non-Gaussianity of the 
lensed sky, but account for the lensing by using accurate 
theoretical lensed CMB power spectra Q . Obtaining pa- 
rameter estimates at Planck sensitivity is then no more 
difficult than with WMAP. 

Fig. |3| shows the posterior parameter constraints from 
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FIG. 3: Parameter constraints from a simple Planck-like sim- 
ulation. Solid lines analyse the lensed sky assuming Gaussian- 
ity with the lensed CMB power spectra, dashed lines are for an 
unlensed sky analysed with the unlensed power spectra, dash- 
dotted lines show the (inconsistent) result from analysing the 
lensed sky using the unlensed power spectra. The bottom 
row shows the dark energy density, age and Hubble constant 
derived parameters. 



this particular realization. The constraints from the 
lensed sky analysed using the lensed Ci are in quite good 
agreement with those one would have obtained if the un- 
lensed sky were observable and was analysed with the 
unlensed C\. So current tools appear to be sufficient to 
extract parameter constraints reliably at up to Planck- 
like sensitivity. Of course it is essential to model the 
lensing consistently: Fig. |3| also demonstrates the wrong 
parameter constraints that are obtained if the lensed sky 
is incorrectly analysed using the unlensed Ci . 

We appear to agree with Ref. [24j that at this sen- 
sitivity lensing does little to change the effectiveness of 
parameter estimation, with the recovered error bars be- 
ing of approximately the same width as when analysing 
an unlensed sky. However we have not tested the accu- 
racy of the recovered error bars, merely showing that the 
posterior peaks are at about the correct parameter val- 
ues consistent with the error bar. The result of Ref. 0] 
would suggest that the error bars are likely to be correct 
to < 5%, but for a precision estimate of the error bars 
a fuller analysis would be required. Increasing the error 
bars should make the lensed and unlensed analysis re- 
sults even more consistent, though there is no reason in 
general why the posterior means should be identical in 
any given realization. 

The effect of non-linear evolution of the lensing poten- 
tial can be safely neglected at Planck sensitivity. How- 
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ever we have made one important linear-theory assump- 
tion, which is that the unlcnscd CMB power spectra can 
be worked out accurately. Given a standard cosmology 
and an ionization history this is straightforward to do, 
however uncertainties in the complicated details of re- 
combination may mean that the ionization history is not 
very well known 5 [2^, |2(| |2tJ ■ Future CMB parameter 
analyses, including Planck, may require a more accu- 
rate calculation of recombination to obtain unlensed (and 
hence lensed) power spectra accurately enough. In ad- 
dition, other second order signatures such as from the 
kinetic Sunyaev-Zel'dovich effect will probably have a 
small but important contribution that also needs to be 
accounted for [H HJ, El| . 



IV. CONCLUSIONS 

We have described a lensed CMB sky simulation code 
that can be used for testing future analysis methods. As 
a simple benchmark we have demonstrated parameter 
estimation at up to Planck sensitivity by using the lensed 
CMB power spectrum on scales I < 2000. This recovers 
the correct parameters, though the error bars may be 
less accurate. More sensitive observations will require a 
fuller analysis of the non-Gaussian likelihood function. 
Higher resolution observations may require a much more 
complicated numerical analysis of non-linear evolution of 
the potential and other non-linear effects |29l |32j . 

Incidentally we have also demonstrated that current 
parameter estimation methodology is sufficiently accu- 
rate for parameter estimation with Planck (under the 
assumption of linearity and assuming that the ionization 
history is known well enough). The simple lensing anal- 
ysis considered here only causes an order unity increase 
in computing time compared to an unlensed analysis, so 
there is no problem accounting for the lensing effect in pa- 
rameter analyses. Future CMB parameter studies must 
account for CMB lensing to obtain correct results. The 
lensed CMB simulation code is publicly available 6 , as 
is the parameter estimation code including support for 
lensed CMB power spectra 7 [T^ . 
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APPENDIX A: HARMONICS AND MAP 
MAKING 

We use spherical polar coordinates, with orthonormal 
basis vectors at any point on the sphere given by eg and 
e^. A complex spin-s quantity can be defined as the com- 
ponents of a rank-|s| tensor in the basis e+ = eg ± ie^. 
For example the spin two polarization [33l | is given in 
terms of the polarization tensor P a \, by 2P = e+ e b + P a {,. 
Further definitions and derivations using our notation 
can be found in the appendix of Ref . [34| . 

A spin s quantity s rj has harmonic components s ai m 
given by 

sV = s ai m s Yi m (Al) 

lm 

where I > \m\,l > \s\. The spin harmonics are defined 
by 

s Y lrn = s \i m e vm<t > = i/ 7~j S ||j 3 s Y lm , (A2) 

and 3 ~' s ' = (— l) s 9 ' s ', Yi m is a standard spin zero har- 
monic and 3 is the spin raising operator (see Ref. |34|). 
The harmonics have the symmetries 

s-Mm — I 1 ) -s x l(-m) 

s Y lm {ir-0,4>) = (-l) l+m _ s Y lm (9,<b). (A3) 

We define gradient (E) and curl (B) harmonic compo- 
nents of a tensor field with the general definition 

\s\Elm = ( — l)^ ^(\s\ a lrn + (—^Y -\s\ a lm) 

i\ s \Bi m = {-l) H \{\ s \ai m - {-l) s -\ s \ai m ) (A4) 

where H £f m = (-l) m | s |£ ;( _ m) , \ s \B^ m 
(~ l) m |s|-S/(-m) an d (— 1) H is a sign convention. 
This definition ensures that gradient fields are always 
pure E. In general a complex spin field s r], with 
s |?7* = -| s |?7, can be expanded as 

\s\V = (~ l) g y^,(\s\Elm + i\ s \Blm)\s\Ylm 

-\s\V = (-l) H+s Y,(\*\ Elm -i\s\Bi m )-\s\Yi m . (A5) 

lm 

For polarization 2-P = e" e b + P a b — Q + iU where Q and 
U are the Stokes parameters measured in the (ee,e^) 
basis. Where Q and U are instead measured with re- 
spect to a left handed set (eg, — e^) as in Refs [j3 Ull 
we have 2P = Q — iU (this basis defines a right handed 
set about the incoming photon direction). In healpix 8 



a http : //cosmocof f ee . inf o/viewtopic .php?t=174 

6 http : //cosmologist . inf o/lenspix 

7 http : //cosmologist . inf o/cosmomc 
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1.2 conventions 2P = Q + iU, and H — 1, whereas in 
the conventions of Refs |3J, |33 H = 0. On large scales 
Cf- 5 = (-l) ff+1 |C, T£; |, and the if = 1 convention is used 
for the output from CAMB and CMBFAST so the large scale 
correlation is positive. To be consistent with these codes 
we use the H — 1 convention for numerical work. 



Harmonic transforms 

For CMB lensing we need to generate maps of the gra- 
dient of the potential, a vector field on the sphere. Vari- 
ous first order results (for example the efficient quadratic 
estimators for the potential 0) also require transforming 
to and from spin one and spin three fields. 

To compute the spin s harmonics one can either iter- 
ate the spin s recursion relation (see Ref. 0), or relate 
the harmonics to the spin zero harmonics. Here we give 
the results for the latter approach, which is followed by 
HEALPIX [10| . Transforms of spin zero fields and spin two 
fields (for polarization) are standard 35] and included in 
HEALPIX 1.2. 

Defining 

= i( 5 A im ±(-l) ro s A K _ m) ) (A6) 
and using some results from Ref. [34[ we have 



same time one can also use a relation like 

I sin 9 s+1 wt = (l- s) (m S W* -lcos9 S W±) 

+ {l + s)a lms Wf l _ l)m . (A9) 

To transform to and from a map of a spin field i s m = 
R + il (where R and / are real) we use the symmetry in 
9 so that 

R(9, 0) = (-if £ isFL\s\E lm + i.F^ [tl B lm ) e m <>> 

hit 
lm 

J2(sFL\s\Ei m - ls F lmlsl B lm )e^ 

lm 

J(tT-0,0) = (-l)H+l+rn+s x 

Y,(° F L\s\Blrn + isF^ A E lm )j™*. 
lm 

(A10) 

For pixelizations where pixels are on lines of constant 
latitude (such as healpix, igloo [sil and glesp 
the <p transform can be performed rapidly using FFTs. 
The remaining computational cost is easily parallelized 
by sending pixels at different 9 to separate processors, 
and Z max « 2000, n p i X « 10 7 maps can be generated in a 
few seconds over about fifty modern processors. 



iW lm = 

= 

2W lm = 
3^+ = 

3^" = 

where 



sin ( 

m 



[a/rnA(;_i) m - I COS 9Xl r , 



sin ( 



2(m 2 - I) 



sin 



1(1-1) 



2 cos 9 

Mm H . 2 a a linX(l~l) r . 



sm 



Otlm\(l-X)m — (i — 1) COs9\lr, 



2m 



sm 



Ki -\)(i- 2) _ 4 2i+m y 3 ) ) ( ,,^ A; „ 

sin 9 



aim \ + 1) + 6- 
m 2 -3Z + 2 „, 



sm 



+12a ; , ( 



4(2 + m 2 )\ 


sm a j 




-l)(i-2)) 


A/m 


cos6»A (; _ 1)m 


m 


sin 2 9 


sin 9 ' 



Oil r, 



I (21 + 1) (F — m 2 ) 



21 - 1 



sinf 



(A7) 



(A8) 



Similar results can be derived for higher spins if desired. 
If several different spin transforms are being done at the 



The deflection field 

The gradient of the potential (a scalar), VV>, has har- 
monic components 



x E lm = (-l) H+1 ^l(l + l)4, lm iB im = (All) 
which follows from 



J dn(e + -VTj>) x Y? m 



- / dn S x Y* m = -y/l(l + l)7p lm (A12) 



= - / dQ 3 if) -iYi™ = (* + l)^m (A13) 



-ia;r, 



where e±- V0 are the spin ±1 components of V0. Thus 
maps of the gradient field can easily be constructed from 
the harmonic components of "0 using results from the last 
section. For spin one field if) = R + il = e + ■ X, we see 
that R and I are simply the eg and components of 
the vector field X. 
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Accurate lensed map making 

To make accurate lensed temperature maps on the full 
sky we use 



f(n)=T(n') 



Im 



T lm Y lm (n') 



(A14) 



Given a set of harmonic coefficients, constructing the 
lensed map is straightforward. Using identities for spheri- 
cal triangles, the lensed temperature at a position is 
given by the unlensed temperature at position (#', <f>+A(j)) 
where 



cos 9' — cos d cos 9 — sin d sin 9 cos a 
sin a sin d 



sin A<j) = 



sin 



(A15) 
(A16) 



and the deflection vector is d = V-0 = dgeg + d^e^ = 
d cos aeg + d sin ae^ . Except near the coordinate singu- 
larities the obvious Euclidean results are a rather accu- 
rate approximation. 

For the polarization the points move the same way, 
with the polarization maintaining the same orientation 
relative to the deflection vector at the two points (ne- 
glecting held rotation, see Ref. [Tl|h However for the 
components of a spin field we have to be careful to ac- 
count for the different direction of the coordinate vectors 
at the two points || . This requires rotating the compo- 
nents of the spherical polar coordinates by 7, the differ- 
ence between the angle made by eg and the connecting 
geodesic at the two points. If a' = a — 7 we have 



tan(a') 



d sin d cot 9 + dg cos d 



(A17) 



and after weak lensing a spin s field s rj becomes 

s ?j(n) = e i ^ 8 r ] (h'). (A18) 

For the spin two polarization field we can avoid inverse 
trigonometric functions by using 

2l7 2{dg + d^A) 2 _ 2i{dg + d 4> A){d 4> -dgA) 
' ~ d 2 (l + A 2 ) d 2 (l + A 2 ) 

(A19) 



where A = tan a' . Except near the poles this is very close 
to unity. 

Note that to get an accurate simulation of the lensing 
B modes at I > 1000 it is necessary to include relatively 
high I (more than just l max + 500 that is accurate for the 
other spectra) . The B lensing signal only becomes useful 
for parameter estimation after Planck. 

To generate the lensed field we can not use FFTs be- 
cause even if n is sampled equally in cf> on rings of con- 
stant 9, the original positions n' will not. An i max w 
2000, ripi x w 10 polarized map can be made in about 
2000 CPU hours, with good trivial parallelization. In 
practice a much faster way to make maps to good ac- 
curacy is to generate an unlensed sky at higher reso- 
lution and just re-map the points. We find that using 
12 x 8192 2 ~ 10 9 pixels for the high resolution map is 
sufficient to get a 12 x 1024 2 -pixel lensed sky with Ci accu- 
rate to 0.5% at I < 2000. The fractional accuracy on the 
polarization is similar, except for the B modes induced by 
lensing which are a more sensitive: 12 x 16384 2 ~ 3 x 10 9 
pixels are required to get Cf B accurate at percent level 
for I < 2000. The B-spectrum is quite sensitive to l max 
and the non-linear power spectrum, so this accuracy level 
is only notional. Note that choosing Z max too low gener- 
ally underestimates the B-mode power, whereas pixeliza- 
tion errors lead to an overestimation, so it is possible to 
get spuriously accurate power spectra without actually 
having computed the lensed sky accurately. 

The approximate point-remapping method is also eas- 
ily parallelized, and by only generating sections of the 
high-resolution map on each cluster node, the total mem- 
ory requirement per node can remain below one gigabyte 
as long as enough nodes are available. The scaling is then 
approximately the same as making maps at the higher 
resolution, and even multi-billion pixel remappings can 
be done in under an hour. For Planck resolution obser- 
vations, lensed simulations can be done in a few min- 
utes with enough processors. For applications where this 
speed is an issue some kind of faster interpolation scheme 
might be useful (avoiding the need to generate a fine map 
at much higher resolution than the base pixelization) . 

The parallelized code incidentally also provides a fast 
method for performing spherical harmonic transforms 
(without lensing) on computer clusters using MPI. 
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